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Abstract. We consider self-similar measures /i with support in the interval 



< a; < 1 which have the analytic functions {e*' 



: 0, 1, 2, ... } span a 



dense subspace in (fi). Depending on the fractal dimension of fi, we identify 
subsets P C No = {0, 1, 2, ... } such that the functions {e„ : n G P} form an 
orthonormal basis for {fi). We also give a higher-dimensional affine con- 
struction leading to self-similar measures fi with support in . It is obtained 
from a given expansive p-hj-u matrix and a finite set of translation vectors, 
and we show that the corresponding (fi) has an orthonormal basis of ex- 
ponentials e'^'^^'^, indexed by vectors A in My, provided certain geometric 
conditions (involving the Ruelle transfer operator) hold for the affine system. 
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1. Introduction: Fractal Measures 



, or in complex form e„ = z" where z — 



We use the notation e„ {x) :— e 
gi27r£c^ a; G R. RecaU that if /i is Lebesgue measure on / = [0, 1] ~ R/Z = T, 
then the functions {e„ : n € No} span the Hardy space H2 of analytic functions on 
T, and {e„ : n g Z} is an orthonormal basis for (I) with normalized Lebesgue 
measure. Functions of the form F (2) = X^i^o ^n-^" with X^J^o ^ °° nis-Y be 

viewed also as defined on U := {z G C : \z\ < 1} and F (z) may be identified with its 
values on the boundary dU = T, represented by z = e*^'^^, a; € R. Writing / (x) :— 



F and defining E, (x) = •£n=o 

kernel formula 



we get the familiar reproducing 



(1.1) 



F{z) = {Ez\f) 



/ E,{x)f{x)dx^y^z''A^, 

•^0 n=0 



\z\ < 1. 



Since we show that, for some fractal measures on / = [0, 1], there are associated 
subsets P C No with orthogonal and total {e„ : n e P} in {n) , it follows that this 
property then carries over to (/i), i.e., the corresponding Hardy space H2 (P, fj) 
is all of (fj,). 

In the present paper, we will consider fractal measures /i supported on compact 
sets (typically totally disconnected) in R'' which arise from iteration algorithms 
that are defined from a given affine system of maps in R''. These maps in turn 
will be defined from an expansive real v-hy-v matrix R, and a given "sparse" set 
of translation vectors in M.'^ . The fractal nature of these constructions arises from 
choosing the number of translation vectors to be strictly smaller than |det R\ where 
R is the given expansive matrix. Our aim will be to understand the possible or- 
thogonal harmonic function systems in (/i), and, in particular, to give conditions 
for (/x) to have an orthonormal harmonic basis. (This continues earlier work 
of the co-authors [IoPe97|.) Many of the ideas going into the general case are 
present already in a relatively simple example in a single dimension, and this will 
be introduced and discussed in Section |^ below. 

While there has been a great amount of recent interest in measures which arise 
from geometric iteration algorithms, see, e.g., |Ban91 , | Ban96[ , |BaGe94|, |BoTa87|, 
Hoffl|, ]LaWa96(| , iLaWa96d| |, ^SS9|, ]Ruc8^ ], and |PoSi95t , there are relatively 
few results on the c orresp onding harmonic analysis of these measures. But recently 
(see, e.g., |BoTa87| and | Hof95 |) diffraction lines in quasicr ystals ha ve lent them- 
selves to such a harmonic analysis. In our earlier paper | JoPe97 |, we noted a 
connection between certain diffraction spectra and the spectral duality considera- 
tions which go into the harmonic analysis of affine self-similarity (for the measures 
fi which we sketched above). 

Wh ile previous harmonic-analysis approaches to self-affine measures (see, e.g., 
Str89t ) have stressed continuous transforms and asymptotic summation methods, 
our present approach is discrete and directly based on Fourier series, for those fi 
which have the orthogonality property under discussion. 

A central tool in our work is a certain double duality: first the usual duality of 
Fourier analysis, corresponding to the dual variables on either side of the spectral 
transform; and secondly a duality which derives from our use of matrix scaling. 
Small scales correspond to compact attractors of fractal Hausdorff dimension, while 
large scales ( "fractals in the large" ) correspond to a discrete set of frequencies (in 
v dimensions), A = (Ai, . . . , Ajy) G R" which label our Fourier basis of orthogonal 
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Figure 1. Support of ^ 

exponentials e\ (x) :— e*^^'^'^ where x is restricted to the ("small scale") fractal. In 
our setup, both scales, small and large, are finitely generated, referring to two given 
finite subsets B and L in W (one on each side of the duality) which are paired in a 
certain unitary matrix U {B, L), defined from the two sets, and to be described in 
( ^^ ) below. The unitary matrix U {B, L) is related to one studied by Hadamard. 
It turns out that not all configurations of sets B, L allow such a unitary pairing, 
and there is a further constraint from the dimension v of the ambient Euclidean 
space. 

We now turn to the theory, following an illustrative example in a single dimen- 
sion, i.e., V = 1. 

2. An Example: Scale Four 



It is known (see, e.g., [JoPe93b|, [IoPe96|) that there is a unique probability 



measure /x on M of compact support such that 

(2.1) [fdt^=U[f(^) dt^i^)+ [f(^ + hdt^i^: 



2 V7 V4y ' J ' \4: 2 

for all continuous /. In fact, the support if of /x is the Cantor set obtained by 
dividing / = [0, 1] into four equal subintervals, and retaining only the first and 
third. (See Figure 1 below.) The Hausdorff dimension of ii is = |^ = i. 

This is a special case of a more general construction in v dimensions {v > 1) 
corresponding to some given real matrix i?, and a finite subset B C M'^. It is 
assumed that 

(2.2) R has eigenvalues all satisfying > 1. 

The subset B is required to satisfy an open-set condition: Introduce 

(2.3) abx ^ R-'^x + b, a; e M". 

It is assumed that there is a nonempty, bounded open set V such that 

(2.4) \JabVcV 

with the union disjoint corresponding to distinct points in B. Our present {at} 
systems (see below for the axioms) are special cases of iterated function systems 



i.f.s.) considered in [Hut81|, see also |Fal86|. There are many interesting more 



general i.f.s., and that context also leads to measures fi which satisfy a general 
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version of the invariance property (2.5), and there is then a corresponding "open- 
set assumption". But for our present ajfine systems, the sphtting property (2.4), 
for some open subset V in R", can be shown in fact to be automatic, see [ JoPe96|. 
In concrete examples, e.g.. Figures |]-^ below, the choice of such an open set V is 
often apparent from the geometry of the afhne maps ab- Wc do not directly use the 
open-set condition ( ^.4| ) in our present proofs, except in the application of Theorem 
3.3 (Remark ^.5[ ) where the interior of a particular simplex is identified for the open 
set V, and the property gets used in ( ^.10 ). For general i.f.s., the property is needed 
for the computation of geometric invariants, such as the Hausdorff dimension, see 
p^alsej , for the particular iteration limits; and, for our present systems, it forces 
the invariant fractal measure /i to be locally translation invariant, see [JoPe97, 
Appendix]. U N — ^ (B), then the corresponding measure /i on M.'^ (depending on 
R and B) has compact support, and satisfies 



(2.5) 



-E 

beB 



f {(Jb{x)) dn{x) 



for all continuous /. (For more details on the "open-set condition" and affinely 
genera te d frac tal me asures, we give the following background references: [ JoPe94 |, 
Str94[ , |Str95[ , and [lPcd97| .) Define, for teR", the Fourier transform 



(2.6) 



AW 



i27Tt-X 



d/i (x) 



with t ■ X ~ Yli=i ^i^ij then get 

(2.7) ^{t)^XB{t)fi{R*-H) 
where 

(2.8) XB(i):=^E^™ 



beB 



and R* is the transposed matrix. For the above example, this amounts to 

1 .wx . /t' 



(2.9) 



A(i) 



;i + e"*) A 



t G 



3. Orthogonal Frequencies and Fractal Hardy Spaces 



Assume that the matrix R in (2.3) has integral entries, and that 

(3.1) RB cZ", Og B, 

but that none of the differences b — h' is in 1/ when b,b' G B are different. Fur- 
thermore, assume that some subset L GlTL^ satisfies € i, # (V) — N {= ^ (B)), 
and the matrix 

(3.2) Hbl := {( 



i2nb-l\ 



is unitary as an iV x complex matrix. 



i.e., Hgj^Hgj^ = In ( * = transposed conjugate). 

In fact, it can be checked that the assumed non-integrality of the differences b — b' 
(when ^ 0) follows from assuming that Hbl is unitary for some L as described. For 
our purposes, the assumptions L Gllf and RB C may actually be weakened as 
follows: 



(3.3) 



(i?"6) • ? e Z, Vn e N, 5 e B, ; e L. 



DENSE ANALYTIC SUBSPACES IN FRACTAL L^-SPACES 

Details on this more general setup will be given in Section |^ below. 
Lemma 3.1. With the assumptions, set 

(3.4) P := {^0 + -R*^i + ■ ■ ■ : U (z L, finite sums} . 

Then the functions {e\ : A € P} are mutually orthogonal in (/z) where 

(3.5) ex (x) := e^^^^'^. 

Proof. Let A = ^ R*^h, ^' — J2 R*^^'i be points in P, and assume A ^ A'. Then 
(ga |eA'),, = J e\ey dfj. 

= J eM^'~^y-d^Iix) 

= A(A'-A) 

= A(/^,-/o + ^* + 



If /q 7^ Zq then xb [Iq ~ ^o) — by (3^). If not, there is a first n such that Z'^ 7^ /„, 
and then 

A (A' - A) = A {R* " (C - z J + i?* (C+i - + • • • ) 

^ XB {I'n ~ L) A {I'n+l ~ L+1 + ■■■) 

= 

since xb (Ki - = 0. □ 



Corollary 3.2. Let /i be the measure on the line M given by and with Haus- 

dorff dimension dn — \ - { We have i? = 4, _B = |0, and L = {0, 1}.) Then 

(3.6) P = {Zo + 4^1 + 42^2 + •••:/, e {0, 1} , finite sums] , 
and {ex : A G P} is an orthonormal subset of {^). 

Proof. Immediate from the lemma. □ 

Lemm a 3. 3. Let the subsets B,L G , and the matrix R be as described before 
Lemma O. Let 

(3.7) Qi(t) ^ lA(i-A)^ teR\ 

Then {ex ■ A G P} is an orthonormal basis for L^ (/i) if and only if Qi = l on . 

Proof. If {eA : A G P} is an orthogonal basis for (/i), the Bessel inequality is an 
identity when applied to et; that is, 

^ = ¥^t\\l = Y.Yex\et)J\ =Y,\fi{t^\)\\ 

A A 

Conversely, if this holds, and if / G (^) {ex : A G P}, then {et \ f) ^ — G for all 
t G K'^, or equivalently / g-^^'^*'^/ (a;) d^i {x) = for aU t <eW . This implies / = 
by Stone- Weierstrass applied to the compact support supp (/i). □ 
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Theorem 3.4. Let H2 {P, 11) he the closed span in I? (/i) 0/ the functions {e'^'^"^ ; 
n = 0,1,4,5,16, 17,20,21,. ..} {i.e., P = {Iq + Ah + A% + ■ ■ ■ : k e {0,1}, finite 
.sums}). Then 

(3.8) ff2(P,/i) 



Remark 3.5. While (/i) is spanned by {2 



e P} (z = e*2^^), and P C No, it 



is of course not the case that 



(3.9) 



z" d/^ = for all n eN. 



(The coefficients do vanish for n E {1, 4, 5, 16, 17, . . . } = P\ {0}, even when n is of 
the form n = 4^ (2Z + 1), fc € No, but not for n in a bigger subset of N.) First of aU, 
we showed in |JoPe9(;| that {z" : n G P} is maximally orthogonal; and secondly, if 
fi is viewed as a measure on T, then the validity of (3.9) would imply (by the F. 
and M. Riesz theorem) absolute continuity of /x with respect to Lcbesgue measure 
A, i.e., ^ & (A), which is clearly false. It is known (e.g., JoPe94| ) that /i is 
purely singular. 



Corollary 3.6. There is a canonical isometric embedding $ of {^) into the 



subspace H2 (z^) + ZH2 (z*) of H2 where H2 (z"*) 
given by 



{/(z4) :/ei?2}; and it 



IS 



(3.10) 



E 

AGP 



c\ex J = C4„z 

neP 



Proof Since EaspI'^^P < 00, and P = [jie{o,i} ^ + 4^', with 4P n (1 + 4F) = 0, 
the representation (3.10) is well defined. Note that $ is everywhere defined on 
(/i) by Theorem 3.4, and the two functions /o (z) = EiiepC4n2*" and /i (z) = 
I]«ep C4n+iz'*" are in (z"*). □ 



Remark 3.7. {Fractal Hardy spaces) An iteration of the argument from the proof 
of the corollary yields, for each n g N, a natural isometric embedding of {p) 
into the subspace of H2 characterized as n increases by: 

H2 (z*") + zH2 (z4") + Z^H2 (z*") + z'H2 (z*") 

+ z'^H2 (z^") + z^'H2 (z4") . . . + Z^i72 (z^") . 

Specifically, let n e N be fixed, and let P„ = {k + Ah H h4"-i;„_i : e {0, 1}}. 

Then the functions in <3>„ (X^ (/i)) (c -ff2) have the following characteristic module 
representation: 

{j2zPf,{z'"):f,eH2}. 
For each n, maps mio this space, and not onto. 



Remark 3.8. {Spectral pairs) Our interest in the problem of finding dense analytic 
subspaces in (/x), for probability measures, grew out of our earlier work on spec- 
tral pairs, see, e.g., [ [JoPe91 |, | JoPe92|, JoPe93a| . Consider subsets and A in M'', 
with Q of finite positive Lebesgue measure, and let (O) be the corresponding 
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Hilbert space from the r2-restricted and normalized j/-dimensional Lebesgue mea- 
sure. Let e\ (x) := g*^'^'*''^ on fl. We say that (il, A) is a spectral pair if {e\ : A e A} 



is an orthonormal basis in (fi). It is known |Fug74] that, for i/ — 2, the case 
when is either the triangle, or the disk, does not admit any sets A such that 
(ri,A) is a spectral pair. On the other hand, our work in |JoPe94 showed that. 



when (f2, A) is a spectral pair in v dimensions, then Q is often "generated" by some 
amount of self-affine structure, as described by a system of afhne transforms ab as 
in (ph. 



Since there is a lack of symmetry of the two sets fl and A in a spectral pair, we 



introduced in |JoPe97| a generalized spectral pair formulation for two measures fi 
and p. The context was locally compact abelian groups: Let G be a locally compact 
abelian group with dual group F. Let /i be a Borel measure on G, and p one on F. 
For / of compact support and continuous, introduce 



{tx)f{x) dp{x) 



where {^,x) denotes the pairing between points ^ in F and a; in G. If / i— > F^f 
extends to an isomorphic isometry (i.e., unitary) of LP' (p) onto {p), then we 
say that (/i, p) is a spectral pair. (The case when p is a restriction of Lebesgue 
measure was studied in |Ped87|.) It is clear how the earlier definition of spectral 
pairs is a special case, even when G is restricted to the additive group M." . But it is 
not immediate that there are examples {p, p) of the new spectral pair type which 
cannot be reduced to the old one. 



Theorem 3.4 shows that this is indeed the case (i.e., that there are examples): 
Let G = M, and let p be the fractal measure in Theorem Let 

P = {Iq + Ah + 42/2 + • • • : e {0, 1} , finite sums} = {0, 1, 4, 5, 16, 17, . . . } , 

and let p — pp he the counting measure of P. Then the conclusion in Theorem 
3.4 may be restated to the effect that {p,pp) is a spectral pair. This is perhaps 
surprising as earlier work on Fourier analysis of fractal measures, see e.g 



Str90c|, 



|Str93 |, and |JoPe95|, suggested a continuity in the Fourier transform, and also the 
presence of asymptotic estimates, rath er than e xact identities. 

It can be shown, as a consequence of [ IoPe97 , Corollary A. 5] that if {p, p) satisfies 
our spectral-pair property for any measure p, then p = pp for some subset P C K", 
i.e., LP (p) has an orthonormal basis of the form {e\ : A 6 P}. The basis for this 
argument is the finiteness of p, when generated from (|2 



Remark 3.9. (The triadic Cantor measure) The significance of the assumptions 
(|^)-(m) on the pair R, B lies in the identity (^^) below, and also in orthog- 
onality. If, for example, we work with the more traditional triadic Cantor set, then 
the results in Lemmas 3.1 and 3.3 no longer are valid. To see this, take i? = 3 and 
-B = {0, |}. Let /i3 denote the corresponding measure on K with support equal to 
the triadic Cantor set (see Figure 0). It is determined by 



/ dp:^ = i 



/ (0 rfM3 {x) + 



/(- + - 
' 3 3 



has d}i — iS-£, and satisfies 

m 3 ' 



/^5(t) 



dp3{x)\, V/eGe(M) 



t e 
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1/3 2/3 



I I I I 

1/3 2/3 1 

Figure 2. Support of ^3 



Choose L = {0, 1} so that (3.2) is vahd; then the subset P3 (which corresponds 
to P = P (L) in Lemma |3.3| ), is 

^3 = 1^ (^0 + 3^1 + 32/2 + •••):?, e {0, 1} , finite sumsj , 

but the corresponding exponentials {e^ : A € P3} are now not mutually orthogonal 
in (/is)- Take for example the two points A = | and A' = | both in the set P3. 
The corresponding exponentials e\ and e\> are both orthogonal to cq, but they are 
not mutually orthogonal, i.e., {e\ \ ey) 7^ 0. In fact, for the /i3-inner product: 



It can further be shown (see | JoPe97 | and Section || below) that there is no subset 



P C K such that, if pp denotes the corresponding counting measure on K, then 
{li3,pp) is a spectral pair in the above more general sense (a fortiori, fractions 



don't provide a basis either). Similarly, it can be checked that the identity ( |4.3| ) in 



Lemma |4T| below fails for this pair /X3, P3, i.e., the triadic one with dn = 



4. Proof of Theorem 3.4: Ruelle's Transfer Operator 



We need only verify the assumption in Lemma But when applied to the 
example in Section ||, the issue becomes showing that J^neP 1/^ (t — n)\^ = 1, t G K, 
where the summation in n is over 

(4.1) P^P{L) = {/o + Ah + 42^2 + •••:/, e {0, 1} , finite sums} . 



More generally, let the sets B,L C K'', and the matrix P, be as in Lemma |3T 
Lemma 4.1. The function 

(4.2) gi(t):=^|A(t-A)|2 



AGP 



{where P = {/q + R*li + ■■■'■ li £ L, finite sums}) satisfies the functional identity 
{t e W) 

(4.3) Q (i) = E {t - l)f Q {R* {t - I)) . 



DENSE ANALYTIC SUBSPACES IN FRACTAL L^-SPACES 



9 



Proof. Let teR". Then 

Q,{t) = J2\f,it-\)f 



xeP 



{t~\)\'\fi{R*-'{t^x))\' 

AGP 

E E - ^ - l'^ ' it-i)-x)\ 
T.\^B {t-i)\'Y.\f,{R*-'{t-i)-x)\' 



leL XeP 



leL 

Remark 4.2. We shall need the following operator C, Ruelle's transfer operator, 
defined on functions Q: 

(4.4) C (Q) it) E {t - Q {R* {t - I)) . 

leL 

We also note 
Lemma 4.3. 



i) The function Qi, defined on by (4.2), has an entire analytic extension to 



C which is of linear exponential growth in the imaginary direction, i.e., with 
C = W +iW, 

(4.5) |Qi(t + is)| <e4™ll^ll2 

{t, s e M.'^) where m is the diameter of the support of fi inside . 



(ii) The operator C in (4.4) extends naturally to the space £,„ of entire analytic 



functions f on C which satisfy (4.5) and f (0) — 1, and maps this space into 
Proof For t e , 

Qi w = E - ^) - ^) ^ E I e*)/. (^-^ I e-*)A. ■ 

AGP ASP 

If the right-hand side is truncated, summing only over the finite subsets in P given 
by the restriction |A| < n, n = 1,2,..., then each of the resulting finite sums 
'^xeP\x\<n' ' ' defines an entire analytic function in t (i.e., t is now allowed to 
range over C, with the convention e* (x) — e'^'^*'^, and t ■ x — Tl!j=i ^j^jj ^ 



Xj G M). These functions are also in L (/i) since the measure /x from (2.5) has 
compact support supp (/i) in W^. Hence 



E ■••|< El---|^EK^A|et)||(e-A|e-,)| 

A 



AGP AGP 
|A|<n |A|<n 

|2\ 2 



A A 

L2(^) ||e-tlL2(^) 
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Note this holds for the P summation, for aU n. Letting n — > cx), we then also get 
the estimate for the unrestricted summation over A G P. This may be formulated 
as a pointwise approximation 

lim g„ (i) = g (<) , <ec^ 



the limit now providing an extension of Q ( • ), initially given only on M.^ by (4.2). 
Suppose the support of /i is contained in {x € : |a;| < m}; then 



|et|lL2(p) 



d/i (x) 



I a; I < rn 



where 



ilmtll = (^|Imi 



2\ 2 



J=l 



t £ 



This means that the functions Qn ( • ) are uniformly bounded in t e , restricted 
to each strip |jlmt|| < constant. We then conclude by the theorems of Montel and 
Vitali (see |Neh75, p. 143]) that the limit Q (i) is entire analytic, t £ C , and serves 
as an analytic extension of Qi from (4.2), given only on M.'^. 

The technique which we used above in estimating the extended Q {t + is) in 
C = R'^ + iW^ also applies, mutatis mutandis, on the other term, to yield: 

2 ^ AT — 1 \ ^ _47rllb|L llsll o 



(4.6) 



by virtue of (3.2). When combined with ( |4.5D , defining and (4.4), the invariance, 
C: £„i £„i follows. We include the details of proof for estimate {il) (using 
unitarity of the matrix in ( [3.2] )): 

leL leL beB 



beB ' 
beB 



i2TTb-{t+is) 



□ 



beB 



Lemma 4.4. Let fj, be a probability measure with compact support on , 
and let P C M.'^ be such that {e\ : A G P} is orthogonal in (/i). Set Qi (t) = 
X^AeP lA ^ -^)l • Then Qi is on ; in fact, it has an entire analytic exten- 
sion. Let A: L^ (/i) L^ (/i) be the orthogonal projection onto H2 {P, /i). Then 

(i) Qi{t) = \\Aet\\l2^^y 

(ii) ijQi it) \t=o = 0, 

(iii) I^Qi (t) \t=o = \\Axj\\'^^ - \\xjf^, where \\xj\\'^^ = Jx'jdfi{x), and 

(iv) ^Qiit)\t^o^i~ir{\\xnl-\\Axnl)- 

Proof. Compute! □ 
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Corollary 4.5. Let {e\ : X G P} be a set of orthogonal frequencies in (p) 
and let Qi (t) — X^AeP lA ~ ^ ^ ■ Assume L is not contained in a hyper- 
plane in . Then the set P is total if and only if 

2n 



(4.7) 



d_ 



Qi(t)|t=o = 0, j = l,...,iy, n^l,2, 



Proof. By Lemmas |3.3| and 4.4, we only need to note that the monomials are dense 
in (/i), by Stone- Weierstrass, and that a vector f E (/i) is in the subspace 
AL^ (/i) if and only if HA/H^^ = ||/||^;. (For more details, see Section ^ below.) □ 

5. The Example Revisited 

Let u = 1, B ~ {O, ^}, L — {0,1}, and R — A. Let /x be the corresponding 
measure {dn (m) = 5)1 and let 

(5.1) P = P(L) = {la + Ah + A% + •••:/, e {0, 1} , finite sums} . 
Lemma 5.1. The function 



gi(i) = ^IA(t-A)r 



AGP 



satisfies the following identity: 



(5.2) 



T:t 



Qi (O=cos^ \ \Q,\ \ +sin^ hr 



■nt 



t - 1 



Let C he the operator defined by the right-hand side in (5.2). On the convex space 
Wo {Q e ([-i,0] ) : Q (0) = 1}, introduce the metric 

(5.3) d(gi,Q2) :=sup{|Q'i(t)-g^(0| [-^0]}. 

Then (Wo ,d) is a complete metric space, and 



(5.4) d(C(Qi),C(Q2))< [\ + ^ 



d{Qi,Q2), VQi,Q2eW^o. 



Pr oof. First note that the functional identity (5.2) in the lemma is a special case 
of (13) in Lemma 41 above. It is clear that CI — 1 where 1 d eno tes the constant 
function. In the general case of Lemma U, that follows from (fj). 

For Qi, Q2 G Wq, we have Qi (i) — Q2 (<) = /g [Q'l (s) — Q'2 (s)) ds, and therefore 
\Qi {t) — Q2 [t)\ < \d{Qi,Q2)- In particular, (Wo,d) is a complete metric space. 
We have C inducing an operator in the convex set Wo, and 



C{Q)'{t) = 



cos 



Q' 



sm 



Q' 



+ - (- siuTTi) Q 



TTt 



+ — (sin TTt) I 



TTt 



t - 1 

4 

t-1 
4 

i - 1 



+ 2 (sinTTi) 



Q' (s) ds. 
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For Qi,Q2 G Wo, the first term is estimated above by j s\rp^^^_i gj \Q'i (t) — Q'2 {t)\ 

since the two affine transformations t 1-^ j and t 1— > leave the interval [— ^,0] 
invariant. The second term is bounded by 

l~ j^JQ'l{s)-Q'2 {S)\ds 

4 

where we estimate |sin7rt| by ^ in the interval t € ["jiO]. Since the integral is 
bounded by j sup^^^r 1 (,i \Q'^ (s) — (s)|, the result follows. □ 



Corollary 5.2. Since the constant function Q = 1 in [—^,0] satisfies C (Q) = Q 



by (5^),, we conclude that Qi f rom (IS 
fixed-point principle. Theorem 



must be constant in 
3.4| now follows by Lemma p. 



0] by the Banach 



6. One Dimension: Even and Odd Scales 

If a subset P C M is orthogonal in the sense that {ca : A G P} is orthogonal in 
(fi), for some fixed probability measure on M, then so is the reflected set —P, 
and the translates t + P when t e K is fixed. In reviewing orthogonal subsets P we 
shall therefore impose the condition G P for simplicity. 

We now examine the following three cases for /i as they illustrate the nature 
of the assumptions made above, as well as the variety of the possible cases: For 
the subset B = {O, |}, we have considered the two scales R = 2 and i? = 4, and 
also R — 3 separately in (6.10) below. The resultin g m ea sure s 112, / i4, an d fi^ are 
determined by their Fourier transforms, as follows, Q, dJ), and ( l6T0| ): 



(6.1) 



00 



1 



n=0 



n 



I = I I e 

ri=l 
TTt 



COS 



e 



n 



COS 



2" 



TTt 

2" 
sin (TTt) 

TTt 



where in the last step we used a familiar infinite product formula for cos(^). 

The infinite product coincides with the integral e 
Lebesgue measure on /. Similarly, 



i27TXt 



dx 



Since ^i2 IS 



(6.2) 



M4 (0 



00 ^ 



n 



n=0 



TTt 

2~4^ 



(The 113 product will be discussed separately.) We picked the set L — {0, 1} (c Z) 
to satisfy condition (3.2) above. Then from (6.1)-(|6^) we get: 



P2 = {h + 2/1 + 2^/2 + •••:/, e {0, 1} , finite sums} = No = {0, 1, 2, . . . } 



and 



P4 = {?o + 4/1 +42/2 + •••:?., e {0,1}, finite sums} = {0, 1, 4, 5, 16, 17, . . . } . 
We get explicit formulas for the zero sets of the transforms, 
(6.3) Z (/I2) = Z\{0} (= the nonzero integers). 



DENSE ANALYTIC SUBSPACES IN FRACTAL L^-SPACES 



13 



and 



(6.4) Z {Jm) = {4" • (1 + 2Z) : n = 0, 1, 2, . . . } 

(= powers of 4 times odd integers) 

where /i2 is Lebesgue measure, supported on / = [0,1], while is the fractal 
measure with dn ~ \ and support on the Cantor set from Figure |l} In the first 
case, it is classical that H2 (P2, ^^2) is the familiar Hardy space in (^2)- That is 
because {^2) may be identified with of the circle, and we are back to classical 
Fourier series, i.e., H2 spanned by 1, z, z^, . . . in (T). The contrast between the 
two cases is made clear from a comparison of the respective functions Q2 and Q^^ 
(derived from ( |3.7|) ) as follows: 

2 ---,2 



Q2{t)= [jT^it^nf^Y. 

neP2 n=0 



sin TT {t ~ n) 



TT {t — n) 



sin {irt) 



E 



1 



- 



where we provide the following interpretation of the last expression (see [ Art64| 
and [Car95|): Let T (t) be the gamma function, with its mcromorphic extension 
resulting from the functional equation 

r(i)r(i-t) 



sin (nt) ' 



recalling 



for positive t. When t > 0, 

d_ 

dt 



T{t) 



T'{t) 



E 



1 



dx 



is well defined, and negative t is from the mcromorphic extension. It follows that 
(6.5) 



^ / N II . II 2 sin (tt^) , , , 



where A2 denotes the projection in onto the Hardy space i?2- (Note that it is 
not immediate from this formula (6.5) that Q2 (t) is entire analytic.) In contrast, 
for Q4 {t) we have 



1; 



i.e., as noted in Theorem 3.4, A4 = /, or, equivalently, the /i4-Hardy space is all of 
(/i4), and embedded as a subspace of H2 via 

H2 (z^) © ZH2 {z^) . 



(See Corollary 3.6 for details.) 



(6.6) 



CiQ) it) 



Q 



In the fi2 case, the operator C, which is analogous to C for /i4 in (5.2), is 

y 

In the fi2 case, the convex Wq is 
(6.7) 



W^o-{QeCi([-l,0]) :Q(0) = 1} 
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It can be checked that C is then not strictly contractive in Wq relative to the metric 

(6.8) d(gi,Q2)= sup \Q'At)-Q'2{t)\. 

te[-i,o] 

and, as noted, C has fixed points in Wq other than the constant function Q = 



1. The contractivity constant can be checked in this case (i.e., (6.6), (6.7), (3.8) 
corresponding to i? = 2) to be ^ + ^ which is > 1 in contrast to the R — 4 case 

in (5^) where the contractivity constant is i + (< l). This illustrates the 
assumption 

(6.9) #(B)<|deti?| 
which will be placed on the general systems from Section |^ below. 



Set z = 



J2TTt 



and g (z) — cos (nt) = 5 + 3-^ + 1^- Then ( |6.6D may be rewritten 



in complex form as 



c{Q){z)^ J2 9HQH, 



which is a special case of the familiar Perron-Frobenius-Ruelle operator [ |Rue94 1 . 
The presence of the above mentioned two independent solutions (in Wq) to C (Q) =^ 
Q is not predicted by Ruelle's theory. The Perron-Frob enius e igenvalue problem 
C (Q) = Q is actually studied in wavelet theory, see, e.g., | Dau92 , Chapter 6]. Since 
{z G T : g {z) = 1} is just a singleton there is onl y on e periodic solution Q in Wq. 
This is consistent with the second solution from (6.5) above being non-periodic 
We also considered B = {O, |} and R = 3, where the measure /is has dn = 
(see Figure 0), and is given by 



In 2 
In 3 



(6.10) 

with 
(6.11) 



n 

n=l 



cos 



2nt 



Z(a^ 



(1 + 2Z) : n = 1,2,3,. 



The contrast to the earlier case is that the sets L which make (3.2) valid, e.g., 
L = {0, |}, are not contained in Z. It can be checked (see Theorem |6.1| below) 
that (/ia) does not contain orthogonal families {e\ : A G P} when P (c M) has 
three, or more, distinct elements. In particular, there are no infinite orthogonal sets 
of exponentials in (ps)- It can be checked, from our formula for Z (jm), that P4 
(= {0, 1, 4, 5, 16, 17, . . . }) is a unique maximal, orthogonal subset for /i4, containing 
0, and contained in Nq. 

We finally mention that ^3^4 :— fi^ * fi^ (convolution) has P4 as a maximal or- 
thogonal set for (^3,4)- (Recall Z (/l3)nZ (jm) — 0, from ( |6.4| ) and ( 6.11 ) above, 
and /i3^ = /I3 • /I4.) But {e^ : A € P4} is not total in (^3,4); or, equivalently, 

Qit) := J2 \f^iit~n)f 

is not constant on M. 

The possibilities described above may be summarized in the following two the- 
orems (v = 1) for measures on M with affine self-similarity. For both results, we 
have R, b determined as follows: R G Z\ {0}, b G M\ {0}; defining ao (x) = R~^x, 
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cTfc (x) — R ^x+b, recall that the corresponding probability measure /i on R depends 
on both R and b, being determined uniquely by 

(6.12) M= ^ 

when \R\ > 1. 

Theorem 6.1. Let R be an odd integer, R =/= ±1, and let fi be the corresponding 
measure. Then any set of fi- orthogonal exponentials contains at most two elements. 



it)=l[-{l + e^^)^e'^l[cos( — ] 



Proof. Recall that 
A 

n=0 

and therefore 

f R^ 

Zifi) = l—i2Z+l):n = 0,l,. 

If 7j, j — 1, 2, 3, are such that the e^. 's are mutually orthogonal in then the 

differences ji - jj (i ^ j) are in Z (fi). Let Ai =71-72, A2 = 72-73, = 71-73, 
and Xj = -^Ij^ {2zj + 1), Zj E Z. Since 

Ai + A2 = Aq, 

we get 

(2zi + 1) + i?"^ (2z2 + 1) = (2zo + 1) , 

which is a contradiction: the left-hand side is even while the right-hand side is 
odd. □ 

Theorem 6.2. Let R be an even 0) integer, \R\ > 4, and let b £ ^Z\{0}. Let 
zq = and set P := zq {Iq + Rli + R?l2 + ■ ■ ■ ■ h € {0, 1} , finite swrns}. Then 
{e\ : A G P} is an orthonormal basis for L^ (/i). 

Proof. The only part of the proof which is not included in the previous discussion 
is the strict contractivity of the operator C: Q ^ C{Q), from (4.4). It specializes 
to 

C (Q) (t) = cos2 (TTbt) Q (J^^ + sin2 (nbt) Q 

Let J denote the closed interval between and Depending on the signs, we 

can have as an endpoint to the left or the right. Set 

Wo = {QeC'{J):Q{Q) = l} 

and 

d(Qi,02) :=sup|Q'i(0-Qi (t)|. 
tej 

Then the argument from the proof of Theorem |3.4| leads to the following strict 
contractivity for Q C (Q) with an explicit formula for the contractivity constant 
(7<1): 

d(c(Qi),c(g2)) <7rf(Qi,Q2) VQi,g2 e Wo, 
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where 



\R\ 



1 

R' 



(The details of the proof of this contractivity property are the same, mutatis mu- 
tandis, as the corresponding ones in Section ^ for the special case R — 4, b = 
Specifically, one of the terms in the estimation of d (C (Qi) , C (Q2)) will be 



TTb\sm{2TTbt)\ 



IQ'i(s)-02 (s)l ds 



for the case when 



< -|, and t e J. 



The sin-term may be estimated by 



sin ( 27r6 



p^pY J — sin ^'^d the integral by |-^| d{Qi,Q2), and the result 

follows. 

Note that 7 — 7 (-R) is a function only of R, i.e., independent of the second 
constant b. Also note that i? 1-^ I7 (i?)| depends only on \R\. In the range of \R\, 
i.e., > 4, I7 (i?)| is decreasing, starting at jgV^+j (< 1) when \R\ = 4. Finally, 
note that the interval J is chosen to be invariant under the two transformations. 



t '-^ and 1 1 



This completes the proof. 

7. Three Dimensions 



□ 



Let B and L be two finite subsets in M'', both containing 0, having the same 
number of elem ents A^, i.e., # (-B) = # (-^) = N; and assume further that the 
matrix Hbl in ( |3.2| ) is unitary. Also suppose L dl}' . A I'-by-z^ matrix i? will be 
given satisfying the following four conditions: 

(i) R = {a^)\^^^, a,j g Z, 

(ii) C Z^ 

(iii) the eigenvalues of R satisfy > 1, 

(iv) N < \detR\. 

Define 
(7.1) P 



Ub [x] 

Pi 

XB it) 



{lo + R*li 
R-^x + b, 
R*x + l, 



R*% + 



€ L, finite sums} , 



and 



1 gi27rfc-t 
6GB 



If /i is the measure (depending on B and R) satisfying (2.5), then 

fL{t) = XB{t)fL{R*-H) , teM.r 

Clearly then 

P=[j {l + R*P), 
leL 

and it follows from the unitarity of Hgi^ and (|)-([]) that this union is disjoint. 
This is equivalent to the assertion that the representation 



A — ^0 + R*li 



R*H2 + 



ik e L) 
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of points in P (L) is unique. Then we proved in Lemmas 4.1 and 4.4 that 

Qit) :=^|A(t-A)r 

AGP 

is entire analytic, and that it satisfies ( [4.3[ ). 

The attractor X (L) of L is the (unique) compact subset X of M'' satisfying 

i?*X IJ (X - ■ 

leL 

To prove that P (L) , the "fractal in the large" , which is constructed from some 
given scaling system {R, B,L) in R'', provides an orthogonal basis of exponentials 
{e\ : A S P {L)} for (/i), using the Ruelle transfer operator, we must exclude 
the special case when L is contained in a hyperplane in K'^. But this condition is 
in fact necessary as we show in Section S below. Moreover, this restriction on L, 



together with the Hadamard matrix property (3^) for the two sets B and L, is 



stringent, and we showed in [ JoPc96 that it can hold only in dimensions — 1, 3, 
and higher. While we can get examples for v — 2 (t he plane) where the axioms, 
including (3_^), are satisfied, we showed in [ JoPe96t that those planar examples 



will "collapse" down into a line: The possible examples in the plane will either be 
essentially one-dimensional, or they will not have H2 {P {L)) equal to (/i), but 
only a proper subspace; see Section ^ for more details. 

The following scholium for M.'^ will serve as a guiding principle for the results 
in Sections || and || below. We will also there discuss two metrics on the testing 
functions Q, both involving the gradient VQ, but one giving a better contractivity 
constant than the other. The proofs, and the evaluation of contractivity constants, 
will be postponed to Sections || and ^. 

Scholium 7.1. The operator 

C(Q) (<) {t-l)fQipi (t)) 

acts boundedly on functions defined on the convex hull of X (L) . 

Proof. Details to follow: see the proof of Theorem 8.3 below, and Section |9[ □ 



The contractivity property of the scholium refers to the metric on the convex 
hull of X {L) applied to the convex set of functions Q satisfying Q (0) = 1, and 
the proof requires all the listed assumptions, including (|i|)-(|v|) for the matrix R. 

But with weaker assumptions, we can still get orthogonality of {e\ : X E P (£)} 
and uniqueness for the L-expansion: 

Lemma 7.2. With the following assumptions: 

(7.2) L C Z", 

(7.3) R{B)cZ'', 

(7.4) HslH*bl - In, 

we have that the representation X = Iq + R*li + R*^l2 + ■ ■ ■ for points in P (L) 
is unique, and the exponentials {e\ : A € P (L)} are mutually orthogonal in (p) 
where fj, is the probability measure on .subject to 



11 i^/ioa, 1. 
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Proof. The crucial step in the proof of both conclusions of the lemma is the obser- 
vation that 

XB {l-l' + R*{p-p')) = 

whenever I and I' are distinct points in L, and p,p' G P (L). First notice that for 
be B: 

eb + R* {p - p')) = eb (/ - /') eb {R* {p - p')) 

^ eb{l- I') em {p - p) 
= (/ - I') . 



We have Rb e W by ([7J), p - p' e , in fact, P {L) C Z", by (|7j), and 

To prove the uniqueness assertion, it is enough to show that, if 

l-l' + R* {p-p') = 

for some 1,1' e L and p,p' e P (L), then I — V . (The uniqueness will then follow 
by induction.) Suppose not, i.e., I ^ I'; then we saw that 

XB {l-l' + R*{p-p'))^0. 

Since xb (0) = 1, we have a contradiction. 

The orthogonality of e\ and ey for distinct A and A' in P (L) follows from the 
representation 

oo 

fAt) = Y[xB {R*-''t) 

11=0 

applied to i = A — A'. The first nonzero term in the L-expansion of A — A' yields a 
factor in the product of the form 

XB {l-l' + R* {p-p')) 

and we noted that this vanishes whenever I ^ V in L. It follows that {ey | e^)^ = 
/t (A — A') =0 which is the desired orthogonality. □ 

Consider subsets B, L d and an integer matrix R with the properties men- 
tioned in Lemma 3.1. We then have dual transforms 

a*i {t) -.^ R*-^ {t~l), teW 

indexed by Z G L. (Compare to the associated transformations { uh : b G B} in 
(2_^).) Assume that {a^ : Z G L} satisfies the "open-set condition" (2^) for some 
open set V* C K"^ which contains in its closure. We now turn to basis properties 
of the exponentials {e^ : A G P} when P = {Iq + R*li + ■ ■ ■ : k e L}. Note that 
Lemmas 3.1 and 5.3 will apply provided we can get 

Q^C{Q) (t) J2 (t - 01' Q (t)) 
contractive in the convex set 

(7.5) Wo := {Q e (W) : Q (0) ^ 1} . 

Note that some V* may possibly be scaled down to get C: Vq Vq strictly con- 
tractive. This will be studied in detail in Sections H and pl 
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Aus dem Paradies, das Cantor uns geschaffen, soli uns niemand 
vertrieben konnen. — D. Hilbert |Hil26, p. 170] 
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Example 7.3. The conditions are satisfied in the following example (|JoPe96 Ex- 
ample 7.4], the Eiffel Tower, see Figure ||): 



(7.6) 



R 



B = 



L = 




The natural candidate for a subset P C such that {e^ : A G P} is an orthonormal 
basis in (fi) is P = {Iq + 2li + 2'^l2 H : k e L, finite sums}. If A S P, the 



three coordinates A = 
space Wo will be 



arc all in Nq. One of the metrics on the corresponding 
d 



(Q1-Q2) 



d(Qi,Q2) := y^sup 

The contractivity proof generalizes that of = 1 in Section Note (in the 
general case 1/ > 1) 



(7.7) 

since 
(7.8) 



leL 



0, 



leL 



Pick a term in ( [7?7| ) which is positive, e.g., 

d 



XB (Or >0; 



then 



E 



ieL\{o} 



(The significance of (7/7) and ( [7.8[ ) will become clear in Lemma 7^ below.) Make 
line integrals connecting to R*^^ (t — I) for each / £ L\{0}, or even curve inte- 
grals (if V* isn't convex). Hence, replacing Qi — Q2 by Q, we find that 



(7.9) 



d 



is estimated by a constant times 
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The property (7.8) is important for the argument, and we include the (easy) 
proof. It is based on the assumption (3.2) placed on the two given finite subsets B 
and L in R". We now calculate from the left-hand side of ( [7.8[ ): 



leL beB 



leL beB b'eB 

beBb'eB leL 

A^^'E E ^* ib-b')Sb,b' 

beB b'GB 

N-'^^et (0) = AT^^iV = 1, 



beB 

which is the claim in ( |7.8| ). 

Remark 7.4. For the Eiffel Tower example we may take 

(7.10) V^={{ti,t2,h):-a<U<0} 

where a e R+ is determined to get the contractive property on Q 
EieL\XBit~l)fQ{a:it)). Set 

3 



C{Q) it) 



ooV' 



and 



Then 



Wo (a) ^{QeC^ {-a < U < 0) : Q (0) ^ 1} 



|-C (Q) (0 - ^ E (t - l)f |-Q {o\ (t)) + second term, 



leL 



and therefore 



d{C{Qi),C{Q2)) < ^d{Qi,Q2) + second term. 



The "second term" comes from evaluating \xb {t ~ 01 when 

A ^ •• . 



(7.11) XB{t-l) 
We use the estimate 



two minus signs 



-y 



-a dt, 



Q 



< 



Q 



and complete the proof of the contractive property of the operator Q i-^ C (Q): 
We get a strict contraction if a S M+ is chosen (so small that) | + ^ < 1, i.e., 
a < ^. Note however that the set ( 7.10| ) will not be invariant under the affine 
dual maps t i— > (t — I), I € L, unless R is expanded by a certain scale, to 
be specified in the following section. We show there that strict contractivity of 
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Q C (Q) may always be obtained by such a "rescaling". (See the details below 
in Proposition |7.6| .) We also give a different metric on the functions Q which yields 
sharper results. Proposition 7.6 implies that the contractivity is indeed strict for the 

matrix R = (oro),r>3, of Example |7.3|. We then get Q = 1 by the Banach fixed 

V r / ' ' 

point principle. Therefore — „ ^ — ;;^Q (t) \t=o = if m = mi + m2 + > 1. 

dt-^ dt^ dt^ ^ 

Hence x^^x^^x^"" e {P (L) , m) and H2 {P [L) , m) = (m) as claimed. 



For use both in the detailed estimates in Example 7.S, and also in Sections g 
and 1^ below, we need the following lemma. Its statement involves some notation: 
II • II2 for the P-norm in v dimension, || • for the corresponding operator norm 
on u x u matrices, i.e., 

||i?||„p :=max{||i?x||2 H^^H^-l}. 

Finally, let | • l^^ y supremum-norm applied to functions on Y, i.e., taking 

supremum over some simplex Y C M". 



Lemma 7.5. Let the setting be as in Scholium 7.1, corresponding to a given ex- 
pansive (see (2.2)) scaling matrix R, and a finite set L of translation vectors in M", 
such that e L. Define L* = L\{0}. 

(i) Then, for every set of coefficients {Qj/g^, Q £ C, such that J2i<£L ~ ^' 
have 

(t) = (P' (^)) = E (^"' - 0) - Q (R-'t)) 

leL leL* 

where pi (t) := R-'^ {t-l). 

(ii) Since 



Q {R-^ {t -l))-Q {R-h) = - [ [Vi) Q) {R-^ {t - si)) ds, 

Jo 

we get the estimate 

\Sc{t)\ < J2 l^'l / \R-H^i)Q{R-Ht~sl)) \ ds, 



leL* 

and, if ^ is the convex hull of the attractor X {i.e., X is the solution to 
X = i + RX), then, for all t G Y, 

l^c (01 <||i?"'IL,E|c'l 11^112 IIIVQII^L.Y- 

(iii) In case R is diagonalizable with a single positive eigenvalue, then Y (= 
convX, i.e., the smallest convex subset in W invariant under all the pi's) 
is equal to the simplex Y' in generated by the points yi := — {R — I) I, 

leL*. 

Proof. In the formula for Sc (t), substitute cq = — X];eL* Then 

Sc (i) = - E (^"'^) + E (^"' " ')) 

leL' leL' 



Y^ci (Q {R-^{t-l))-Q{R-h)) 



leL' 
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as claimed. The difference terms inside the sum may then be estimated using the 
auxihary functions Lpi (s) Q {R^^ {t — si)), < s < 1, where 

^(p, {s) = -R-'iVi)Q{R-'it~sl)), 

and noting that, for each s G [0,1], I E L, and for t E Y, the points R^^ {t — si) 
arc all in Y. 

To prove the asserted invariance of the simplex Y', let Y denote the convex 
hull (in R'') of the attractor X of the afRne system {pi : I E L} (i.e., the unique X 
(C R") satisfying X = L + RX). Note that, for each I E L, yi = - {R- I)^^ I is 
the unique fixed point of pi, including the case I = 0, where po (0) = 0. 

For n G N, and {li,l2, ■ ■ ■ ,ln) G L x L x ■ ■ ■ x L , let y {h, . . . , In) be the fixed 

n times 

point of pi^o pi^o ■ ■ ■ o pi^ , and let X„ be the set of all the fixed points y {h, . . . ,ln) 
as {li, . . . , In) varies over X ^ Then it is known (see |Hut81, Theorem 1]) that 
X — lJ^gpjX„, and therefore 



(7.12) Y = convX = |J convX„ 



where conv denotes the convex hull. Since po (0) = 0, convXi is the simplex Y' 
generated by {yi : I E L*} where yi — — {R ~ 1)^ I, i.e., 

Y' = convXi = [Y, siVi ■Si>Q, H ^ l} • 



Using now the contractivity property of the affine maps {pi : I E L} implied by (2.2), 
we infer by induction that, if R is diagonalizablc with a single positive eigenvalue, 
then 

(7.13) X„ C Y' for all nEN, 



and, therefore, by (7.12), that Y C Y'. A further volume consideration shows that 
in fact Y Y'. The invariance of the simplex follows, as claimed. 

To prove ( 7.13D , proceed by induction, checking first that X2 C convXi when R 

Cr ■ ■ ■ o\ 
r • • • \ 
: : • . : for some r E N. Every y E X2 has the representation 
6 6 • • • r / 

y = {l - R^)^^ {lo + Rli) for some ^o, h e L. It follows that y = {R + I)^^ xq + 
R{R + I)^ xi where Xi := — (i? — /)~ h E Xi for i — 0,1. So when the scahng 
matrix R has the stated form, it follows that y = 7^2:0 + if^^i £ convXi. Hence 
X2 C convXi. Now suppose (7.13) has been proved for n, and let y E X„-|_i. Then 

y= {I - _R"+i)~^ {I0 + Rh^ h i?"^,i) for some k E L, and it follows that there 

arc z E X„ {z = {I - R")-^{lo + Rh -\ h R'^-^ln-i)), and a; S Xi, such that 

y = /ifl^'+i ^ + i+R+^..+R^ ^- Clearly y^rTTnT + i+r+...+r^ = 1' so that, when R has 
the stated form, it follows from the inductive hypothesis that y E conv (X„ U Xi) C 
Y', concluding the proof. □ 

Note that Y is not in general contained in conv(Xi), even if R is in diagonal 
form. For examples when conv (X) is of a more complex nature, see [ BrJo96b| 
Sections 9-10 and the graphics]. 
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The first application is to the contractivity of the C-operator in Example 7.3 
Let the R in (^^) be variable; i.e., R— ( [j g ^ ^ ; rather than the r = 2 version (7.6) 
and assume r G N, r > 2. Then, for each r, we have a corresponding simplex (as in 
the lemma). It is given by parameters > 0, i = 1, 2, 3, such that si + S2 + S3 < 1. 
Specifically, 

(7.14) 

Y(r)= {t= eM3:ii = ^,<2 = ^,i3 = ^,s.>0Vz,|:s,<l}. 

(See Figure |[) Recall that Y := Y (r) is (by the lemma) then invariant under the 
four affine transformations, 

t^— (leL), 

r 

and also e Y, whenever t e Y and s e [0, 1]. 

Let S and i C be as in ([7.6|) of Example |7.3|. For r G N, let 



Pr (L) = + rli + r^^2 + ■ ■ ■ : li £ L, finite sums} , 

a(/)(t) = ElxB(i-OlV(^), /eC-(M3), 

leL \ ^ / 



and 



}rit):^ \l^it-X)\\ te 

\ePr{L) 



p3 



We then have the following conclusions concerning this example. 

(i) If r is even, then ( |3.3| ) holds, an d w e then have the eigenvalue equation 



Cr {Qr) — Qr Satisfied (by Lemma [4.l[ ) 



(ii) If r > 3, then Cr is contractive in the suprcmum-gradient norm (Proposition 
|7.6| below). Hence, if r > 4 and even, then the exponentials 

Er (L) {ex-.Xe Pr (L)} 

form an orthonormal basis in {^Ji-r)- 

(iii) But (J)-(|i|) do not carry over to the case when r in N is odd. For example, it 
can be checked directly that C3 (Q3) ^ Q3, and that (L) is not orthogonal 
in L'^ifis). (The points (0,0,0) and (4,4,0) are both in P3 (L), but not 
orthogonal in (^^3). Specifically, /ij (4, 4, 0) ^ fi3 (|, |, O) 7^ by a direct 
check.) 



Proposition 7.6. In the Eiffel Tower example {Example 7.3), the operator Q 1— > 
Cr (Q) is contractive when r > 3 on the space of functions Q € C°° (R"^) such that 
Q (0) = 0, and relative to the norm 



\Q\\i,Yi3) --^ sup ||VQ(i)|| 
teY(3) 



where \\VQ\\ := E -=1 4« 
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t3 




Figure 4. The Simplex Y (r) and Its Vertices 



Proof. The first term in Cr (Q) of Example 7.3 is ( |7.9D , and we will restrict, for 
simplicity, to functions Q G (R"^) such that Q (0) ~ 0. For the derivative terms 



JrlXB (<-OI'of (l7j|),wehave 



E 



llmxs (t-OI 



;eL* fe=i 



with the sign-convention (and vanishing of one of the three terms in the sum) from 
(7.11) above. 
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Hence 
(7.15) 



E 



< — 2^ |sin(7rtfe)| < — sm 

k=l 



r-1 



when t £ Y(r). (If r = 2, we bound |sin(7rife)| by 1.) The last estimate is from 
evaluation of the maximum of t he X] fc=i -term when t is restricted to the simplex 



Y (r) from Figure ^; see also ( [7.14| ). (The maximum is attained at one of the 
vertices.) 

Returning to the estimate on the sum ( |7.9| ), we have (using Lemma 7.5): 

't - sV 



(7.16) 



|(|7l)|<-sin 



ds 



where estimate (7.15) was used. But when t is restricted to Y (r), then G Y (r), 
and the term under the integral may in that case be estimated by a supremum over 



Y (r) . The corresponding supremum of 

but a more specific terminology would be | 

( [7.16| ), we get, for i e Y (r), 

TT 



Q 



VQIIi 



will be denoted simply ||V(5||, 
oo Y(r)- Continuing estimate 



< — sm 
- 2 



Since V0I3 (Y(r)) 



3(r-l)^ 



r- 1 
we then get 



ivol3(Y(r))||VQ|l 



< 



2r(r~l 

The other term in the estimate of 
t e Y (r) , is bounded above by ^ 



r- 1 



IVQIl 



(Q) it) 



IS 



all' 



which, for 



,Y(r) 



The identity (7.8) was used a 



second time for that. The same estimates work also for the other partial deriva- 



tives 



dt 



■Cr (Q) (t) , and summing the estimates on X]j=i 



mj^r iQ) 



we get 



|va(g)|| < -livgii 



2(r-l)^ 




Summary. With respect to this gradient-norm, Cr has 7 ■ ^1 + 2{r^i)^ ('^)) 
as operator bound. 

This estimate does not get us contractivity of Q 1-^ (Q) (on functions Q G C°°, 
restricted to Y (r), and satisfying Q (0) = 0), when r = 2; but, when r = 3, we 

! .53 < 1 as an upper bound on the contractivity 



get the number ^ (^l 



2-(3-l)" 



constant. 



□ 



Remark 7.7. Note that, rather than using the simplex Y (r) in the estimates for the 
operator C'r of Example 7.3, we could alternatively have used simply the enveloping 
unit cube B positioned as follows: —1 < tj <0,j = 1,2,3, in R^. (See Figure ^.) 
But its volume is one, as opposed to V0I3 (Y (r)) — j^jr^-jja, which gives the present 
much better contractivity constant for Cr. 
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8. A Scaling Property for Matrix Dilation 

Let the system (R,B,L) be as described, i.e., i? is a v-hy-v integral matrix 
which is strictly expansive. The subsets B and L in R"^ will be assumed to satisfy 



(^)-(|3J), specifically R{B) C Z'^ and L C Z'^. If R is scaled by some r e N, 
then the new system {rR,B,L) satisfies the same conditions. We introduce the 
transformations 

(8.1) arM{x) ■■= {rRy^ x + b, xeR", 

and the corresponding measure /i^ (of course also depending on B) . The fixed-point 
property of fir is 

where N = ^ (B) (= ^ (L)), and we recall that fir — ij, for r — 1. The correspond- 
ing formula for the Fourier transform /I^ is 

it) = XB (t) 1m- {{rRy^ tj , 

where, as before, 

fees 



The transformation corresponding to (|4.3| ) is 

(8.3) a (0) (t) = ^ \XB [t - 01' Q (^"'^* {t - 0) ■ 

Example 8.1. We have considered the following two measures, and ^2, in one 

dimension, arising from B = {O, L = {0, 1}, and R = 2. Then rR = 2 or 4, for 
r — 1 and 2, and we showed that fii was Lebesgue measure on [0, 1] while 112 was 
the measure with dn {112) = \ from Theorem We considered the eigenvalue 
problem Cr {Q) — Q, and showed that Q = 1 is the only solution normalized to 
Q (0) = 1 when r ~ 2, but not when r = 1. 

We will now show that this is general, also for the examples in higher dimen- 
sion. We show that the eigenvalue problem Cr {Q) — Q always will have a unique 
(normalized) solution provided only that r is taken large enough. 

The proof is based on two technical points: 

Lemma 8.2. For the directional derivative Va, a G R"^; of \xb (^)I'j we have 

Va \XB {t)\^ = "27riV-2 a-{b-h') sin [2ti {b - b') ■ t] , t e . 

b,b'eB 

Proof. An elementary computation. Recall the notation Va := a • V = ^iW' 
relative to the canonical coordinate system in . □ 

With the assumptions from above on the system {R,B,L) in M"^, we have the 
following estimate for the contractivity constant 7 of 

C (Q) it) = \XB {t -l)\'Q {R* {t - 0) • 
leL 



28 



PALLE E. T. JORGENSEN AND STEEN PEDERSEN 



Let Y denote the convex hull of the attractor X (L) of the maps 
(8.4) pi{t) := R*~^ {t-l), teR", 

indexed by Z e L. 

We replace the norm || • y from Section ^ by the following alternative norm: 



(8.5) 



IQIIi,Y:=ll|VQl2lli = ^|VQ|2dm, 



where m denotes Lebesgue measure. 

Theorem 8.3 (L^-Theorem). Let (R, B, L) he a system in W satisfying (9.1)- 
( |j.3| )^ e i. Let C be the operator given by (4.4), and let Y denote the convex hull 
of the attractor "X. p. Let \\Q\\i y given by ( ^.5[ ), and let 



/3 := 27rdiam(B) max ||sin(27r(6 - 6')( ' -0)llr 

b.b'eB 

leL 



Let \\T\\^p be the operator norm, and 
Schmidt norm for a i> x i> matrix T . Then we have 

(8.6) IICQIIi Y < \AetR\ 



1/2 



the Hilbert- 



{l-N-')p\\R-\^raa^\l\, + N\\R-\^ 



IQIIi. 



Y ' 



for any -function Q such that Q{0) = 0. If further^ C\{Y — I) is a set of Lebesgue 
measure zero for any I g L with I ^ 0, then we have the sharper estimate 



•7) IICQII, Y < Ideti?| 



IQIIi, 



Y ' 



(I- N-^) B\\R-^\\ maxm„ + 

V / ' W M op ' '.^ M Whs 

for any -function Q such that Q(0) = 0. 

Remark 8.4. Note that this result is pointless in one dimension (y — 1), for then 

|deti?| II 11^^ = 1, 

while this product gets small in higher dimensions. In our | JoPe96 | setup, the L^- 
Theorem can only be used if |deti?| ||-R~^||^5 < 1 is possible for some invertible 
integer matrices R. Actually this may work better if we take the generalized setup 

as — > oo. 



3) seriously, and, e.g., take R = i/j., since then |deti?| ||i? ^| 



hs 



However, it is not yet clear if there are any nontrivial examples allowing us to use 
only the L^-Theorem. 

Proof. We have 



\CQh 



leL 



2lll 



= I |E ■ - 01') Q°pi + Y.\^bC -i)\'^{Qo Pi) 

leL leL 

As before we estimate the two terms in this sum, starting with the first term. Note 
that Lemma 7.5 implies 
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for all t. Setting L* := L\{Q}, it follows that 

E \XB{t - l)f) Q o pi{t) = 5^ (V \XB{t - 01') {Q o pi{t) - Q o po{t)) 

E (viXB(t-Ol') f\R*-'l)-{VQ{R*-H-sR*-H)) ds. 



leL 

for all t. So we have 



< 



< 



E -01') , f\R*-H)-{yQ{R*-' ■ -sR*-H)) 

E -Ol')U^*-'^l2 f\^Q{R*-' ■ -sR*-H)\^^ 



ds 



<max V|xb(- 
i^Lj mi 



maxli?* ^l\ 



E t |V(3(i?*"^ • -sR*-H)\ 

leL' Jo 



ds 



We continue the estimate by estimating the last factor in this product. 

ds 



E \VQ{R*-H-sR*-H)\ 
leL' Jo 



/ E / \VQ {R*-H- sR*-^l)\^ dsdt 
Jy (g^. Jo 

/ E / \VQ{R*~^t-sR*-^l)\^dtds. 

Jo ,^r.Jy 



Making the substitution u = R* — sR* = Ps,i{t) we have dt = |deti?| du, so 
that 



E f\yQ{R-h-sR-H)\^ds = fj2 [ 



\VQ{u)\2 \detR\ duds 



< 



|VQ(u)|2|deti?| duds 



<(7V-l)|deti?| IIIVQI2II1. 

The second to last estimate used the fact that ps^i{Y) C Y, which in turn is a 
consequence of the fact that (Ps,;(i))o<s<i the line connecting pi{t) and po(i), 
the invariance p{Y c Y, and the convexity of Y. So we have shown 



V|xs(- -01' m^^\R*-H\^{N-l)\detR\ \\\VQ\^\\^. 

2 00 lEL* 



< max 
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Estimating the LR* term, 



o Pi 



< (1 - 7V-i)/3 11^^ maxilla |deti?| \\Q\\,^^ . 



dt 



leL 

The other estimate is 

\£\xb{ - ~l)\'W{Qopi% = [ Y.\^B{t-l)f\ViQopi)it)\ 

I ; 

<E / \^ iQ o Pi) {t)\, dt 

I 

<E / \i^Q)opi{t)\,\\R-'\\^^dt 

= |deti?|E/ \iyQ)iu)\,\\R-\^du 

< \det R\N\\\VQ\ J ^\\R-X^. 

The last inequahty used p;Y C Y; and it can be improved if Y n (Y — ^) is a set 
of Lebesgue measure zero for every I E L* . In that case, we have 



\Y.\XBi- -l)f\'^iQ°Pi)\2\l < \detR\\\\VQ\J^\\R-' 



I hs 



It is now easy to complete the proof. 



□ 



Remark 8.5. The following additional arguments may be used in a sequence of 
steps (to be described below) leading to the elimination of the factor |det(i?)| in 
the two estimates ( ^.6[ ) and ( |8.7| ) in the conclusion of Theorem 8.3: After estimating 
II V(5||ii(-Y) by a constant times 



E / l(vg)op,| dt 



for Q E C°° {Ey) such tha t Q (0) = 0, we may change variables in the J., 



terms in the summation (p.SI ). Let Y have the non-overlap property (2^), i.e., 
suppose that the Lebesgue measure of the possible overlap in the union 



(■■■)dt 



(8.9) r(Y):=Up,(Y) 

leL 

is zero. For the sum in (|8.8| ), we then have 

V/ \\7Q\o pidt^\dei{R)\ I \\/Q\dt 

and 

/ |VQ| dt < sup |VQ|vol,(T(Y)) 
Jt(y) t(y) 



(8.10) 



"(Y) 

TV vol. (Y) 
= sup VQ |. ■ 
T(Y) |det(i?)| 
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Then combining estimates ( |8.8D and ( 8.10| ), we get 

Y I l(VQ) opi\dt< TV vol, (Y) sup |VQ| . 
77l T(Y) 



When combined in turn with the estimate from Theorem 8.2, we get 
live (Q)llii(Y) < N sup |Vg| (R) vol, (Y) , 

T(Y) 

where F, (R) is an expression (see Theorem 8.3 and Section ^ which depends on 
the operator norm, and the Hilbert-Schmidt norm, applied to R^^. In the case of 
Example O ("The Eiffel Tower"), [f o r 1 = const, -r-^ 

' ' L V r / . 

When the same method is applied to VC" (Q), it gives a new tool for showing 
lim„^oo C" (Q) = when Q (0) = 0. Since 



(8.11) 



X = f|T" (Y)^ 



and X has Lebesgue measure zero when N < |deti?|, the limits can be estimated 
further. 

Recall that the r"-terms in the intersection ( 8.11 ) are defined from iteration of 

T"(Y)= |Jp^,o...op,„(Y), 

(ii,...,i„)eX?i 



(B.c), i.e 



and also that we are finished after showing C" (Qo) for all Qo e C°° (M") such 
that Qo (0) = 0, because the fixed point C (Qi) — Qi {Qi = 1 Qo) will then be 
unique subject to Qi (0) = 1. 



Remark 8.6. The constant (3 in ( ^.6| ) may be easily calculated for the three-dimen- 
sional example in Figure ||, i.e.. Example 7.3. We have (3 — 27rdiam(i3) = Stt-^ = 

7r-\/2. Calculating det (R) in this example, it follows from the formula for (i?) that 
the number + l) F3 (R) (< 1) is a contractive upper bound for the constant 



r 



0; 



7 when i? = I r ) , and r G N is taken sufficiently large. See formulas (9.15) and 



( j.l6| ) in Section |^ for details. We expect r — 2 will suffice. (Hence we get strict 
contractivity in this new metric as claimed.) An application of Lemma 3.3 then 
shows that (/x^) of Figure ^ has an orthonormal basis {ex ■ A € Pr} when Pr is 
given in (7.1). 



Corollary 8.7. When the contractivity constant is calculated for the system 
{rR, B, L) as described above, then limr^oo 7r = 0, and it follows that Cr is strictly 
contractive when r increases. As a consequence, if L is not contained in a hyper- 
plane in My , then the exponentials {e\ : A G P,. (i)} will form an orthonormal basis 
for (nr) where Pr (L) = {Iq + rR*li + ■ ■ ■ : li ^ L, finite sums}. 



Remark 8.8. The formula for the contractivity constant 7^ contains a factor r~^, 
but we also note that the convex hull Y^ of {L) {— the attractor of {p^^i : I G L}) 
decreases as a function of r. (See the next section for details.) 
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9. Estimates and Duality of Scales: Fractals in the Large 

In this section, we provide a context in W (the dimension v arbitrary) which 
utilizes two types of duality: first the traditional time/frequency duality of harmonic 
analysis, and secondly the duality between the iteration limits when the scale (given 
by a relative integral v-hy-v matrix R) gets arbitrarily small, versus when it gets 
arbitrarily large. 

The axioms below are such that the scaling matrix R is only required to be 
integral with respect to the two given (finite) translation sets B and L; see (3.3), 
or (9.2) below. In this general setup, we give explicit estimates on the contractiv- 
ity constant, making it clear how it depends on the given data (i?, B, L) and the 
dimension v. The estimates are gradient-supremum-norm ones. 

Consider a triple {R, B, L) such that R is an expansive v y. v matrix with real 
entries and B and L are subsets of R'' satisfying: 

(9.1) N:=#B^#L; 

(9.2) R"-b-le Z, for any neN, beB, I eL; 

(9.3) Hbx := N-'^^^ (^'^'''''OheB leL ^ unitary N xN matrix. 

We introduce two dynamical systems 

crfc(x) := R^^x + b; 
Ti{x) := R*x + I 

and the corresponding inverse functions 

ujb{x) := R{x - b); 
pi{x) ■.= R*-\x-l) 



where x , b G B, and leL. We showed in | JoPe94| that one may think of {ab) 
and (r/) as dual dynamical systems. It is known that there are unique non-empty 
compact sets Xo- and such that 

(9.4) X,= \Jat (X,) = B + R-^X, 

beB 

and 

(9.5) Xp=[Jpi{Xp) = R*-\Xp-L); 

leL 

and it is easy to see that 



= {Y^R-^bk : bk e b} and Xp ^ {-^.R* ~% ■ h e L^^ 

fe=0 fe=l 



By uniqueness, these sets are th e closur es of the respective orbits of under the 
maps (cTb), respectively (pi). In [ JoPe9(| , we argued that it is natural to think of 
Xp as a dual of the attractor X^. Let /i = /io- be the unique probability measure 
on such that 

(9.6) p^N-'Y.t'°''b'- 

beB 
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Then this measure has as its support. Similarly, Xp is the support of the unique 
probability measure /i' = /ip satisfying 

(9.7) //=iV-i^A^'op-i. 

leL 

The corresponding expansive orbits 

n 

(9.8) C = Xr:= {^R* ''h-ne N, h G l} 

fe=0 
n 

(9.9) B = X^ := [-Y,R''bk-neN,bkeB'j 

k=l 

of under r and lu are also important. It should now be clear that the pairs {at, v) 
and {pijLOb) share many properties, which however we do not always make explicit 
in the following. Notice that the set C has a recursive property similar to ( |9.6| ): 

(9.10) C = L + RC. 

We discussed uniqueness of the corresponding decompositions of elements in C 
above. We also investigated when the functions {e^ ; A G £} form an orthonormal 
basis for L^{p). One may equally well ask when {e^ : A G S} is an orthonormal 
basis for L^(/x'), but that is an essentially equivalent question. It turns out that 
the measure fi' is important for our investigation. We will see that the set B 
has properties that imply that translates of Xo- tile M.^ , when N = |det_R|. The 
assumptions on (i?, B, L) insure that the functions {e\ : A G £} are orthogonal, so 
we really only have to investigate totality. 



The following result, combined with [LaWa96d, Theorems 1.1 and 1.2], shows 
that, if = |deti?|, then X^. and Xp are translation tiles. We must also check 
that C is uniformly discrete, i.e., that there is an e > such that ||A — A'||2 > £ for 
all A and A', A ^ A', in C. That follows from our assumptions on (i?, _B, L), as can 
be verified by an induction argument. 

Lemma 9.1. Each point X E C has a unique representation A = '^R*''lk with 
Ik G L. Similarly each point in B has a unique representation of the form 'Y^R^b^ 
with bk £ B. 

Proof. Each point A G £ has a representation of the form X = I + R*p with I G L 
and p e C. Suppose I - I' + R* {p - p') =0 where 1,1' G L and p,p' G C. We will 
show that I = I' , the desired result, then follows by induction. If G -B, then 

l = eb{l-l' + R*{p-p')) 

^et{l~l')et{R*{p-p')) 

= eb(l - I )eRb{p-p ) 

= eb{l-l'), 

where the last equality used (|9.2l) . If / 7^ /' then we have X^bes — I') = by 



(3.2). But this contradicts (|9.1lD because this sum would be N if (9.11) holds 



hence 1 = 1' as we had to show. □ 



We will show that C has the basis property by showing that C from (4.4) acts 



contractively on a set of functions including Qi of (p.7[) and 1 (where 1 denotes 



the constant function) and thereby showing X^app I/* ( ' ~ ~ That is, we 
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want to use the uniqueness part of the Banach fixed-point theorem on a suitably 
chosen metric space. It is natural to consider C as an operator on a set of smooth 
functions Q satisfying Q(0) = 1. However, the appropriate choice of a metric is not 
so obvious. The present discussion is based on gradient-supremum-norm estimates. 

Let H2{C) denote the subspace of L'^ilj) spanned by the orthonormal set {ex : 
A G £}. Any et, t e C, is in L^{fJ,), so H2{C) is a subspace of L'^{p)- We will 
show that H2{C) — L'^{n) for certain systems {R,B,L) satisfying (|9.lD-(p^). Let 
A denote the orthogonal projection of onto H2{C), that is 

Af := E {ex I /), ex 
xec 

for any / G L'^il^)- Our next result is a useful criterion for when H2{C) = 

Lemma 9.2. Let dj :— If £ L, then H2{C) = L^ifJ-) if and only if 

d"' ■ ■■d'J'^QiiO) = for all > with ni + ■ ■ ■ + > 1, where Qi (t) = 
EAeplA(<-A)|'. 

Proof. Restating the definition of Qi, we have Qi{t) = ||Aet||^ for all t. If H2{C) = 
L'^in), then Aet — et for all t, hence Qi{t) = 1 for all t E M.'^, in particular, 
9"^ • • • 9""Qi(0) = for all nj > with ni + ■ ■ ■ + > 1. Conversely, suppose 
9"^ • • • d"''Qi{0) = for all nj > with rii + • • • + ri^ > 1. By Lemma 4.3 we have 



a power series expansion Qi{z) — J2ne% ^nz"; and, by our assumption, a„ — 
unless n = 0. Hence Qi is a constant function. Since € L, we have Qi{t) ~ 1 
for all t E M.'^ . It follows that et S H2{C,) for all t G M"^, so an application of the 
Stone- Weierstrass theorem implies that H2(C) contains all continuous functions on 
Xct, hence H2{C), being a dense subspace of must be equal to L^{n) as we 

wanted to show. □ 



Let Y denote the convex hull of the attractor given by ( |9.5[ ), let HQHq^ '■— 
suPygY \Qiy)\, and let 

(9.12) IIQIIy^-IIIvqUIoo 

1 /2 

where jzjj := l^jl"^^ is the usual Hilbert norm on C. We begin by finding 

an explicit formula for the operator norm of C acting on a suitable set of smooth 
functions. 



Theorem 9.3. Let {R,B,L) be a system in W satisfying (9.1)-(9.3), G L. Let 



C he the oper ator given by (4.4), let Y denote the convex hull of the attractor Xp 
given by ( |9.5| ), and let \\Q\\y be given by ( 9.12 ). If L spans , and if there exists 



7 < 1 so that ||CQ||y < tIIQIIy f'^'^ '^^^ Q in a set of -functions containing 
1 - Qi, where Qi (t) = Y^xeP - A)l^ H2{C) = L^{fi). 

Proof. Since CI = 1 and CQi = Qi we have 

||1-0i||y = ||C(1-0i)||y<7I11-Qi|Iy 

so Qi{y) = 1 for all y € Y. Since Y is a convex set containing 0, and Y is not 
contained in any proper subspace of M.'^, it follows that 9"^ • • • 9," '^Qi (0) = for all 



Uj > with ni + • • • + rtjy > 1. Hence an application of Lemma 9^ completes the 



proof. □ 
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Remark 9.4. In the statement and proof of Theorem |9.3| , Y could be replaced by 
any set containing Xp. But Y will be chosen so as to produce optimal estimates. 

The restricting assumption in Theorem |9.3| (and also in the corresponding results 
in Sections |^ and |7|) that the set L C W not be contained in a hyperplane is 
necessary. In poPe96 , Example 7.3], we considered a triple {R,B,L) in where 
v = 2, N ~ 3, and where the conditions in Theorem 9.3 all hold, except that L, 
in this |JoPe9(:| example, is contained in a line in the plane. Specifically, take 



6 ^ 
6 > 



B 



R 

- , . . 

measure on M satisfying 



{(o).(S)'(l)}'^ = {(o)-±^i}where 



Properties 



3) are clearly satisfied for the example. Let /i be the unique probability 



A (*) = XB {t) A 



t e 



where 



XB {t) = -{l + e"'^+e"'^) 



Let P (L), Qi, and C be defined (as usual): 

P (L) = {lo + 6li + e'^h H ■■he L, finite sums} , 

Qi(t):= IA(^-A)|^ 



xeP{L) 



and 



C(/) {t) = J2\xB it-l)ff 



t-l 



t e 



From Lemma 4.1, we get the eigenvalue identity C {Qi) — Qi satisfied, and con- 
clude that the exponentials E (P (L)) = {e\ : A S P (L)| are orthogonal in (/z). 
Moreover, the argument from the proof of Proposition 7^ gets us contractivity for 
C, i.e., 

(9.13) ||VC(/)||y<7||V/||y 

for some constant 7 < 1, and for all / £ C°° (M^) such that / (0) = 0; where 



|v/||. 



sup (^|9,/(t)|), 



and Y is the L-attractor. But Y is a finite line in K by Lemma 7.5; in fact. 



Y = {( (2 ) : ^2 = —ti, \ti\ < j^} ■ From the uniqueness of the Banach fixed point, 
we get that Qi = 1 on Y, but not on M^. Indeed, if we had Qi = 1 on M^, then it 
would follow from Lemma |9^ that H2 {P (L)) — (fi), i.e., that E (P (L)) is an 
orthonormal basis for (/i). But it is not, as the function fj_ (x) = e'(^i+^2) is in 
L'^{^), but not in H2{P{L)), since ( } ) -L P{L). That would make f± constant 
on every line of slope 45 degrees and therefore constant on X (B). But we can see 
that it is not from the formula for /i and its support, the attractor X(i?). This 
attractor has Hausdorff dimension see Figure |^, and the formula for X (B), 



(9.14) 



X(P) 



i=0 



1 



h e B 
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X2 




Figure 5. The Cut and Projection Fractal X (i?) 



From (9.13), and Lemma 4.3, we conclude that Qi has the form Qi{ti,t2) = 
F {ti +t2) for some entire analytic function F of one variable. If E {P (L)) were an 
orthonormal basis for (/i), it would follow similarly that every f £ LP' (/i) had 
the representation / {xi,X2) = 3/ {xi — X2) for some function gf of one variable, 
which we noted is not the case: that would make f±_ {xi, X2) = e' ^'^'-"'"^''^ constant 
on X {B) which it is not, as can be checked by inspection of ( 9.14 ). 

The argument shows that functions / € H2 (P {L}) are characterized by the 
following representation: let 



r,, ■.= {X-h:XeP{L)}. 



Then f{xi,X2) = J^^eTi^ c (0 {xi - X2), where J2^er,^ |c(OI^ < 



Inspection of ( 9.14 ) yields the following representation of points x — {^D in 
X(B): 



1 °° / e \ 1 
1=0 ^ ^ 
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where £i,?/i G {0, 1}. Let 7r± (x) xi ± X2, and consider the hnes l± (h) given by 
xi ± X2 = h, respectively, for /i G M, \h\ < 1. Expanding h in base 6, 



1 °° 1 



2 ^ 6 

we have, for the intersections X± (h) := X (B) n l± {h), that X-|_ (h) 7^ if and 
only if hi € {0, 1}. In that case 

X+ (h) = e X (B) : TT+ (x) = /i} 



e {0, 1} 



i=0 

Hence both H := {h E M. : X+ (h) 7^ 0}, and the X+ {h)'s, are fractals (embedded 
in one dimension). The respective Hausdorff dimensions are (H) = and 
Dh (X+ (h)) < 1^, with Dh (X+ (/i)) = i^, generically 

On the other hand, the other line-intersections are at most singletons, i.e., 
X_ (h) = {cc G X (B) : 7r_ (x) — h} is empty, or a singleton. We have X_ (/i) ^ 
if and only if the digits (hi) are restricted to hi e {0, ±1}; now a set of Hausdorff 
dimension . The singleton intesection is then 



X {h) - |x - (77,, eO e X(B) : 77, = if /i, = 0, 

and rji = 1, Si = — if hi e {±1} 

Summary. The Hardy subspace H2 {P (L)) in (/i), /or the X (S) fractal in R^, is 
canonically associated with the following one- dimensional pair: (tt- (X(i3)) jF/j^), 
where 7r_(X(i3)) is f/ie above mentioned fractal in {of the same Hausdorff 
dimension |^), and Fj^ is the associated spectrum. 

O ur n ext result shows that one can find a bound for the constant 7 in Theo- 
rem 9.3, the bound depending only on the data (R,B,L) in a very explicit way. 

Recall, is the operator norm, and :— k=i l^jM ) is the Hilbert- 

Schmidt norm for a // x // matrix T. 



Theorem 9.5. Let {R,B,L) be a system in satisfying (9.1)-( |9.3| ), £ L. Let 

C be the operator given by (4.4), and let Y denote the convex hull of the attractor 
Xp given by (9.5). //HQUy given by (9.12) and 

(9.15) 7 := {N-lfN-'p\\R-\^nm^\l\, + \\Rr\^, 
where 

/3 27rdiam(B) max ||sin(27r(6 - 5')( • -0)IL' 

b.b'eB 

leL 

1/2 

ll^llop ^he operator norm of R, and \\T\\f^^ := (^X]}'fe=i l^j-fel^) « x 

matrix T , then we have 

(9.16) IICQIIy <7||QIIy> 
for any -function Q such that Q(0) = 0. 
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Proof. Begin by observing that 



b.b' 



7V-2^cos(27r(6-6')(t-/)) 



b,b' 



SO 



V {\xBit - 01') I =N-'J2 1^ - ^'l2 \MMb - h'){t - 0)1 

^ b.b' 



Since the N terms with 6 = 6' are 0, it foUows that 



(9.17) 



^ / 2 oo 



for each / e L. Since p/Xp C Xp and pi is affine we have pfY C Y. If < (3(i) < 1, 
then clearly we have the following estimate: 

(9.18) IICQIlY-lllVCQyU 



Yy - - if) Qopi + Y. ■ - if V(Q o 



< ||El^(ix^(- -or 



Qo pi 



first term 



second term 



Let us begin by estimating the first term in this last sum. Using Lemma 7.5 we 
have 

J2^{\xBit~l)f)Q{pi{t)) 
I 

= -T.^{\^B{t~l)\^) r {R*-H)-WQ{R*-'{t-sl)) ds. 
leL' •'^ 

Since (i?*^^(t — sO)o<s<i is the hne connecting poit) and pi{t) it follows from ( |9.17 ) 
that the right hand side is dominated by 

(9.19) ||E|v(ixs(- -or 



1^0 



Qo Pi 



<{N-1){1~ N-') p ||i?-i||^^max \\Q\\^ 



The (N — 1) factor is from an application of Lemma 7.5, In the estimates above 
we used the invariance p; Y C Y and the convexity of Y. 



To bound the second term in ( 9.1q ), we begin by noticing that Lemma 4.1 and 
(^J) lead to 



(9.20) 



|El^s(- -Ol'|V(Qop,)l2 <max|||V(QopOl2l 
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If Cj denotes the jth column in R* ^, then 

d 



Qo pi = ((VQ) o pi) ■ Cj 



so that 



<mQ)opi)\, 



1/2 



1/2 



•J l2 



and it follows that 

||Eixb(' 



< IIQIIyII^" 

Ol'|V(QopOl2 



I /is 



IQII 



Y ' 



hence the proof is complete. □ 
Corollary 9.6. Let {R,B,L) satisfy ^)-{^, and, for r EE, let 

n 
k=0 

and let fir be the probability measure solving 

fi^ ^ N-'^ ^ Pr ° crr,b 

beB 

where Cr^bix) (tR) ^ x + b. Then {e\ : A G Cr} is an orthonormal basis for 
Lp'{pr) for r sufficiently large. 



Proof. This is an immediate consequence of Lemma 3.3 and Theorem |9.5| since 
-1 



\rR\\-; ^ and 



{rR)-' 



as r ^ cxD. 



□ 



hs 



In one dimension, with R an integer, and N = 2, the general situation is sum- 
marized in the following strengthening of Theorem |6.2| above. 



Corollary 9.7. Suppose v ^ I, N ^ 2, B ^ {0,6}, with b e M\{0}, R is an 
integer with \R\ > 2, and p is given by ( |6.12| ). If R is odd, then L'^{p) does not 
have a basis of exponentials for any b G M\{0}. If R is even, and \R\ > 4, then 
L^(^) has a basis of exponentials for all b e M\{0}. If R > 4 and even, then L may 
be chosen in such a way that L^{p) is identified as the boundary values of analytic 
functions in the unit disc. 

Proof. The argument for odd R is in the proof of T heo re m |6.l| . 

Suppose R is even, and let L := {0, i^jy }. Then (|t])-(|]3^ are satisfied, and we 
can easily find the convex hull, Y, of the attractor Xp. (See Lemma 7.5(ii|).) In 



fact, if i? > 0, then Y is the closed interval with endpoints 



(|26|(ii-l)) 



and 0; while 



if i? < 0, then Y is the closed inteval with endpoints (|2fc|(_R^-i)) 719777717^7: 



(\2bR\{R:^-l))- 



Using the notation from Theorem 3.5, we have |ji?| 
and 

/3-27r|&|||sin2^fe.|U. 



op 



\Rl \\R- 



I hs 



R- 
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Prom the characterizations of Y, it follows that 

TT 



/3 = 27r|6| 



sm 



R-1 



and 



P = 27r|6| sin ■ 



i?2 _ 1 ' 

Using N = 2, and maxjgi 6 jljj = j^, we conclude that 



2R 



sm ■ 



R- 1 



and 



7 



2|i?| 



sin 



1 



1 



if i? > 0, 



if < 0. 



if > 0, 



if i? < 0. 



Thus 7 < 1 if I J?| > 4, but not iiR = -2, 2. 



□ 



If i? = 2, then {e^ : A e Z} is an orthonormal basis for L'^in), and clearly 

{eA : A G £} cannot be an orthonormal basis for L'^ip.) for any choice oi L = {0, 1} 
with ^ ^ 0; in fact, depending on the sign of I, the set £ will consist solely of 
non- negative real numbers, or solely of non-positive real numbers. We do not know 
what happens for general b ^ 0. 

If i? = —2, then we do not know if Lp'ilj) has a basis of exponentials for some 
h e M\{0} . 
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